Evidence for a J/ippp Pauli Strong Coupling ? 
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The couplings of charmonia and charmonium hybrids (generically 9) to pp are of great interest 
in view of future plans to study these states using an antiproton storage ring at GSI. These low 
to moderate energy typp couplings are not well understood theoretically, and currently must be 
determined from experiment. In this letter we note that the two independent Dirac (7 M ) and 
Pauli (o>i/) pp couplings of the J/tp and if)' can be constrained by the angular distribution of 
e + e~ — > {J/tp, tp' ) — > pp on resonance. A comparison of our theoretical results to recent unpolarized 
data allows estimates of the pp couplings; in the better determined J/tp case the data is inconsistent 
with a pure Dirac (7^) coupling, and can be explained by the presence of a cr MJ/ term. This Pauli 
coupling may significantly affect the cross section of the PANDA process pp — » ir°J/ip near threshold. 
There is a phase ambiguity that makes it impossible to uniquely determine the magnitudes and 
relative phase of the Dirac and Pauli couplings from the unpolarized angular distributions alone; we 
show in detail how this can be resolved through a study of the polarized reactions. 

PACS numbers: 11.80.-m, 13.66.Bc, 13.88.+e, 14.40.Gx 



I. INTRODUCTION 

Charmonium is usually studied experimentally 
through e + e~ annihilation or hadronic production, 
notably in pp annihilation. The pp annihilation process 
was employed by the fixed target experiments E760 and 
E835 at Fermilab, which despite small production cross 
sections succeeded in giving very accurate results for 
the masses and total widths of the narrow charmonium 
states J/ip, ip' , Xi an d X2- A future experimental 
program of charmonium and charmonium hybrid pro- 
duction using pp annihilation that is planned by the 
PANDA collaboration Q at GSI is one of the principal 
motivations for this study. 

Obviously the strengths and detailed forms of the cou- 
plings of charmonium states to pp are crucial questions 
for any experimental program that uses pp annihilation 
to study charmonium; see for example the predictions 
for the associated production processes pp — » 7r 0, F in 
Refs.@,ii]. (We use 'F to denote a generic charmo- 
nium or charmonium hybrid state, and tp if the state has 
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Unfortunately these low to moderate energy produc- 
tion reactions involve obscure and presumably rather 
complicated QCD processes, so for thepresent they are 
best inferred from experiment. In Ref.[4| we carried out 
this exercise by using the measured pp partial widths to 
estimate the coupling constants of the J/tp, tp' , Vc, v'a Xo 
and xi to pp, assuming that the simplest Dirac couplings 
were dominant. These >Fpp couplings were then used in 
a PCAC-like model to give numerical predictions for sev- 
eral associated charmonium production cross sections of 
the type pp — ► 7r 0, F. 

In this paper we generalize these results for the J/tp 
and vb' by relaxing the assumption of 7 M dominance of 
the vbpp vertex. We assume a ippp vertex with both Dirac 
(7^) and Pauli couplings, and derive the differen- 

tial and total cross sections for e + e~ — ► i\> — ► pp given 
this more general vertex. Both unpolarized and polarized 
processes are treated. 

A comparison of our theoretical unpolarized angular 
distributions to recent experimental J/tp results allows 
estimates of both the Dirac and Pauli J /ippp couplings. 
There is a phase ambiguity that precludes a precise de- 
termination of the (complex) ratio of the Pauli and Dirac 
J ' /ippp couplings from the unpolarized data; we shall see 
that an importance interference effect between the Pauli 
and Dirac terms leads to a strong dependence of the un- 
polarized pp — > ir°J/ip cross section near threshold on 
the currently unknown phase between these terms. 

Determining these couplings is evidently quite impor- 
tant for PANDA, and can be accomplished through stud- 



2 



ies of the polarized process e + e~ — > J / '?/> — > pp. The 
angular distribution of the unpolarized, self-analyzing 
process e + e~ — > J/-0 — > AA may also provide comple- 
mentary information regarding the closely related J/ipAA 
vertex. Both of these processes should be accessible at 
the upgraded BES-III facility. 



II. UNPOLARIZED CROSS SECTION 

The Feynman diagram used to model this process is 
shown in FigJT] We assume a vertex for the coupling of 
a generic 1 vector charmonium state ip to pp of the 
form 



2m 



(1) 



In this paper m and M are the proton and charmonium 
mass, r is the charmonium total width, and we assume 
massless initial leptons. Following DIS conventions, q v 
is the four momentum transfer from the nucleon to the 
electron; thus in our reaction e + e~ — *• pp in the cm. 
frame, we have q = (— t/s, 0). The couplings g and K are 
actually momentum dependent form factors, but since 
we only access them very close to the kinematic point 



q A = M 2 in the reactions e + e 
will treat them as constants. 



pp, we 
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FIG. 1: The Feynman diagram assumed in this model of the 
generic reaction e + e~ — > ip — > pp. 

The unpolarized differential and total cross sections for 
e + e~ — > ip — > pp may be expressed succinctly in terms of 
the strong ippp Sachs form factors Ge — + Ks/Am 2 ) 
and Qm = g(l + k). Both Qe and Qm are complex above 
pp threshold, in part because phases are induced by pp 
rescattering. If we assume that the lowest-order Feyn- 
man diagram of FigfTJis dominant, the phase of g itself 
is irrelevant, so here we take g to be real and positive, n 
however has a nontrivial phase. We express this by intro- 
ducing a Sachs form factor ratio, with magnitude p > 
and phase x\ 



Qe/Qm = pe l 



(2) 



The corresponding relation between the Pauli coupling 
constant n and this Sachs form factor ratio is 



k e 



%4>K 



pe 



1 



(M 2 /4m 2 - pe**) 



(3) 



where we have assumed that we are on a narrow reso- 
nance, so we can replace s by M 2 . 
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We will first consider the unpolarized process e 
ip — > pp, and establish what the differential and total 
cross sections imply regarding the ippp vertex. The un- 
polarized total cross section predicted by FigQ]is 



47TQ 2 M 4 (1 - Am 2 / s) 1 / 2 
3f 2 ~ [(s - M 2 ) 2 + T 2 M 2 } 



(2m 2 \g E \ 2 + s \g M \ 2 ) 



(4) 

(We use angle brackets to denote a polarization averaged 
quantity.) Exactly on resonance (at s = M 2 ) this can be 
expressed in terms of the ip partial widths 



and 
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(l-4m 2 /A/ 2 ) 1 / 2 
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which gives the familiar result 



(2m 2 \g E \ 2 +M 2 \G 



M 



\s=M 2 



12-k 
~M 2 



B, +f 



Bpp . 



(5) 



(6) 



(7) 



and tp — » pp 



Here B e + e - and B p p are the ip — > e" 
branching fractions. 

Since the (unpolarized) pp width and total cross sec- 
tion on resonance involve only the single linear combina- 
tion (2m 2 |t/£| 2 + M 2 |C/m| 2 ), separating these two strong 
form factors requires additional information, such as the 
angular distribution. The unpolarized e + e~ — > ip — > pp 
differential cross section in the cm. frame is given by 



M 4 (1 - 4m 2 /s) 1 / 2 



4/2 



s 2 [(s-M 2 ) 2 +r 2 M 2 ] 



^m 2 \g E \ 2 {\- p 2 ) + s \g M \ 2 {i + p 2 ) 



(8) 



where p — cos{9 c . m .). This angular distribution is often 
expressed as 1 + ap 2 , where 



i-{4m 2 / s )\g E /g M \ 
i + {4m 2 / s )\g E /g 



(9) 



M 



Inspection of Eqs. (j8|9[ ) shows that one can determine the 
magnitude p — \Qe/Qm\ of t ne Sachs form factor ratio 
from the unpolarized differential cross section, but that 
the phase x of /5m is unconstrained. 

The undetermined phase \ implies an unavoidable am- 
biguity in determining the magnitude and phase of the 
Dirac and Pauli ippp couplings g and k from the unpo- 
larized e + e~ — > (J '/ip, ip 1 ) — > pp angular distribution. We 
will discuss this ambiguity in the next section. 
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III. COMPARISON WITH EXPERIMENT 
A. Summary of the data 

Experimental values of a have been reported by several 
collaborations. The results for the J/ip are 



1.45 ±0.56, 

1.7 ± 1.7, 

0.61 ±0.23, 

0.56 ±0.14, 

0.62 ±0.11, 

0.676 ± 0.036 ±0.042, 



and for the ip' 



0.67 ±0.15 ±0.04, 
0.85 ±0.24 ±0.04, 



MarkI 
DASP 
Markll 
Marklll 
DM2 9 

bes m 



E835 
BES [ 



(10) 



(11) 



For our comparison with experiment we use the sta- 
tistically most accurate measurement for each charmo- 
nium state, and combine the errors in quadrature. This 
gives experimental estimates for a of 0.676 ± 0.055 and 
0.67 ± 0.155 for the J/ip and ip' respectively. 



B. Testing the pure Dirac hypothesis 

We first examine these experimental numbers using the 
"null hypothesis" of no Pauli term, k = 0, in which case 
a = (1 - r)/(l + r), where r = 4m 2 /M 2 . This k = 
formula was previously given by Claudson, Glashow and 
Wise [lH and by Carimalo [lj]; the value of a under 
various theoretical assumptions has been discussed by 
these references and by Brodsky and LePage [15], who 
predicted a = 1. Fig[2] shows these two experimental 
values together with the pure Dirac (7^) formula for a. 
The ip' case is evidently consistent with a Dirac (7^) ip'pp 



coupling at present accuracy, but the better determined 
J / tp angular distribution is inconsistent with a pure Dirac 
J/ippp coupling at the 4a level. 

The discrepancy evident in Fig [2] may imply the pres- 
ence of a Pauli term (k 7^ 0) in the J/ip pp vertex. Inspec- 
tion of our result for a in the general case (Eqj9j) shows 
that one can certainly accommodate this discrepancy by 
introducing a Pauli term. 



C. Determining p= \Qe/Gm\ from a 

The dependence of the predicted a on p at the J/ip 
mass (from EqJ5]) is shown in FigJ3] The experimental 
value a — 0.676 ± 0.055 (shown) is consistent with the 
Sachs form factor magnitude ratio of 



P = \Qe/Qm \ = 0.726 ± 0.074. 



(12) 
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FIG. 3: The experimental value of the unpolarized e + e~ — > 
J/ip — > pp angular coefficient, a — 0.676 ± 0.055 (shaded), 
and the resulting Sachs J/ippp strong form factor magnitude 
ratio p= \Ge/Gm\ (Eq[l2j. 
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FIG. 2: The coefficient a observed in the unpolarized e + e~ — > 
(J/ip, ip') — > pp angular distributions, together with the the- 
oretical result a — (1 — r)/(l + r) predicted by a pure Dirac 
(7n) V'PP coupling. 



In terms of p and \ this completes our discussion: 
Given the unpolarized angular distribution, one obtains 
a result for p = \Qe/Qm\ from Eq[9l but the phase x of 
Qe/Qm is undetermined. However one may ask the more 
fundamental question of what values of the Dirac and 
Pauli coupling constants g and k in EqQ] are consistent 
with a given experimental unpolarized angular distribu- 
tion. 



D. Determining k, 

First we consider the experimentally allowed values of 
k. The unpolarized angular distribution provides us with 
a range of values of p (Eq fT2|) , but x is unconstrained; we 
may combine this information through EqJ5]to determine 
the locus of allowed (complex) values of k. This is shown 
in Fig[U 
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FIG. 4: The locus of complex k (the J/ippp Pauli coupling) al- 
lowed by the experimental constraint p — 0.726±0.074, taken 
from the unpolarized differential cross section for e + e~ — > 
J/ip ->pp. 
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FIG. 5: The value of the overall Jjippp vertex strength g 
implied by the experimental Tj/^^ p p and p as a function of 
the unknown J/tppp Sachs phase \. 



For x — 0, Eq[2]implies that k is real and negative, and 
takes on the smallest allowed magnitude. As we increase 
X from 0, the allowed k values proceed clockwise, since k 
initially acquires a negative imaginary part. The extreme 
values of k on the real axis in FigfJ]are for x = 0, tt, and 
are 



-0.137 ±0.032, x = 
-0.500 ±0.011, x = tt- 



(13) 



E. Determining g 

Next we consider the determination of the overall 
J 'ftp — ► pp vertex strength g. Since the differential and 
total cross sections for e + e~ — > J/tp — * pp only involve 
g through the ratio g/f^, we must introduce additional 
experimental data to constrain g. The partial width for 
J/ip — > pp is especially convenient in this regard, since it 
only involves the strong J /ippp vertex, and thus depends 
only on g and k (and kinematic factors). This partial 
width was given in terms of the strong Sachs form fac- 
tors in Eq|5J as a function of g and k it is 



i 



- + 33?(k). 



1 + =- 



This generalizes the k = result given in Eq.27 of 
Ref. [H to a nonzero Pauli coupling. Using the PDG 
values 16] of Tj^ = 93.4 ± 2.1 keV and Bj/^^ pp = 
(2.17 ± 0.07) • 10~ 3 , EqfTHimples a range of values of the 
overall vertex strength g for each value of the (unknown) 
phase x- This is shown in Fig [5] There is a range of un- 
certainty in g at each x ( n °t shown in the figure) , due to 
the experimental errors in Tj/^p, Bj/^_, p p and p, which 
is at most ~ ±5%. 

Note that g is bounded by the limits at x = and tt, 
for which g w 2.0 • 10 -3 and f» 3.4 • 10 -3 respectively. 
The allowed values of g are somewhat larger than our 



previous estimate of g = (1.62 ±0.03) • 10 -3 [4j assuming 
only a Dirac J/ippp coupling, as a result of destructive 
interference between the Pauli and Dirac terms. 



IV. EFFECT ON a{pp n° J/ip) 

The effect of a J /ippp Pauli term on the pp — ► tt J/ip 
cross section may be of considerable interest for the 
PANDA project, since one might use this as a "cali- 
bration" reaction for associated charmonium production, 
and the Pauli term may be numerically important. Al- 
though we have carried out this calculation with the ver- 
tex of EqQ] for general masses, the full result is rather 
complicated; here for illustration we discuss the much 
simpler massless pion limit. 

For a massless pion the ratio of the unpolarized cross 
section (a(pp — * ir°J/ip)) with a Pauli term to the pure 
Dirac result (7^ only, denoted by D) is 




FIG. 6: The dependence of the unpolarized, near-threshold 
cross section (a(pp — > -K°J/ip)} on the (complex) Pauli cou- 
pling k = |«|e**» (from EqlTSJ. 



where j3 = y/l — Am 2 / s is the velocity of the annihilat- 
ing p and p in the cm. frame. The limit of this cross 
section ratio at threshold is shown in FigJS]for a range of 
complex k. 

Evidently there is destructive interference for a n with 
a dominant negative real part, as is suggested by the un- 
polarized data. For the value K — —0.50 (the larger solu- 
tion in Eq [T3")) there is roughly a factor of two suppression 
in the cross section over the prediction for a pure Dirac 
coupling. The suppression however depends strongly on 
the phase of n, and for imaginary k has become a moder- 
ate enhancement. Thus, the near-threshold cross section 
for pp — > 7T J/ip is quite sensitive to the strength and 
phase of the Pauli coupling; it will therefore be impor- 
tant for PANDA to have an accurate estimate of this 
quantity. In the next section we will show how both the 
magnitude and phase of k can be determined in polarized 
e + e~ — ► J/tf) — > pp scattering, and may be accessible at 
BES. 



V. POLARIZATION OBSERVABLES 

The relative phase x °f the J/ippp Sachs strong form 
factors Ge and Gm may be determined experimentally 
through a study of the polarized process e + e~ — > J/ip — > 
pp. As each of the external particles in this reaction has 
two possible helicity states, there are 16 helicity ampli- 
tudes in total. All the helicity amplitudes to the final pp 
helicity states |p(±)p(±)) are proportional to Ge, and all 
to |p(±)p(=p)) are proportional to Qm- In the unpolar- 
ized case these are squared and summed, which leads to 
a cross section proportional to a weighted sum of \Ge\ 2 
and |<?Af| 2 - As we stressed earlier, this implies that the 
phase x of Ge / Gm is not determined by the unpolarized 
data. 

To show how x can be measured in polarized scat- 
tering, it is useful to introduce the polarization observ- 
ables discussed by Paschke and Quinn ^7[- These are 
angular asymmetries that arise when the polarizations 
of particles are aligned or anti-aligned along particular 
directions. For example, for our reaction e + e _ — > pp, 
Q(0, 0, z, 0) is the difference of two angular distributions, 
(da/dQ)pf — (da/dQ,) p ^. Here we will use x and y for the 
two transverse axes and z for the longitudinal axis (see 
Fig|7]). x and z vary with the particle, and y is chosen to 
be common to all. An entry of signifies an unpolarized 
particle. Since there are four possible arguments for each 
particle, 0,x,y and z, there are 4 4 = 256 polarization 




FIG. 7: Axes used to define the polarization observables. 

observables for this process. Of course there is consider- 
able redundancy, since they are all determined by the 16 
helicity amplitudes. The constraints of parity and charge 
conjugation reduce this set to 6 independent helicity am- 
plitudes, and for massless leptons (as we assume here) 
this is further reduced to 3 independent nonzero helicity 
amplitudes. 

We introduce the normalized polarization observables 
Qe e+ e B _e f e p = <3(e e + , £ e~ , e p, <h>) I Q(0,0,0,0), where 
Q(0, 0,0, 0) is the unpolarized differential cross section. 
The (nonzero) polarization observables for this process 
satisfy the relations 



(a) 


Qoooo = 


Qxxyy — 


Qyyyy = 


— QzzOQ = 1, 


(b) 


QoOyO = 


QxxOy = 


QyyOy = 


QzzOy — 




— Qoooy = 


QxxyO — 


~QyyyO — 


— QzzyO, 


(c) 


QxxOO = 


QyyOO = 


QoOyy = 


— Qzzyy, 


(d) 


QzOxO = 


QzOOx — 


Qyxyz — 


Qyxzy — 




— QozxQ = 


— QozOx = 


Qxyyz — 


Qxyzy 7 


(e) 


QzOzO 


QozOz — 


Qxyxy — 


Qyxyx — 




G.i i = 


— QzOQz = 


Qyxxy — 


Qxyyx 7 


(/) 


Qxxxx — 


Qyyxx — 


Qzzzz — 


— Qaozz, 


(9) 


QoOxx — 


Qzzxx 


Qxxzz — 


— Qyyzzi 


(h) 


QoQxz = 


Qxxzx — 


Qyyzx = 


Qzzzx — 




Qxxxz — 


— Qyyxz — 


Qzzxz — 


— QoQzx, 


(i) 


QxyxO = 


QxyOx = 


QzOyz — 


QzOzy — 




QyxxO 


QyxOx — 


— Qozyz — 


— Qozzy 



(16) 



Explicit expressions for these observables are given in 
Table D 

The results in Table [J suggest how we may determine 
X experimentally. Inspection of the table shows that only 
four of the independent polarization observables depend 
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TABLE I: Nonzero inequivalent polarization observables in 
e+e~ -> J/ip -> pp. The function F is 4 - 2(1 - p 2 ) sin 2 0. 



Pol. Observable 



Qoooo 

QoOyO 
QxxOO 
QzOxO 
QzOzO 
QoOxz 
QxyxQ 
Qxxxx 
Qzzxx 



Result 



1 

4p sin x sin 9 cos 9/F 
2(1 -p 2 ) sin 2 0/F 
4p cos x sin #/f 
4 cos 9/F 
— 4p cos x sin 9 cos 0/F 
— 4p sin x sin 0/F 



[4 - 2(1 - 
-2(14 



i 0]/F 
•'6»/F 



would suggest plausible J /ippp couplings. This approach 
has some experimental advantages; as the A and A decays 
are self-analyzing, no rescattering of the final baryons is 
required to determine their polarization. In addition no 
beam polarization is required, since it suffices to measure 
the (odd-p) polarization observables Qoo y o and Qqo xz - 
One may also measure the even-p observables Qooxx and 
Qoozz as a cross-check of the result for p. 

Finally, we note in passing that it may also be possible 
to measure the appropriate polarization observables in 
the time-reversed reaction pp — > J/ip — ► e + e~. 



VI. SUMMARY AND CONCLUSIONS 



on x; two are proportional to sinx and two to cosx- As- 
suming that one knows p with sufficient accuracy from 
the unpolarized data, one may then determine x unam- 
biguously by extracting sinx and cosx from the mea- 
surement of two of these polarization observables. 

The determination of sinx is the most straightforward, 
since it only requires the detection of a single final po- 
larized particle (for example the proton, through Qooyo)- 
If k is close to real, which corresponds to x ~ or w tt, 
this observable may be relatively small. The other polar- 
ization observables that are proportional to sin \ involve 
asymmetries with either one or three particles polarized; 
these are given in relations (b) and (i) of EafTol 

Determining cosx involves measuring double or 
quadruple polarization observables, which are given in 
relations (d) and (h) in Eq[T|)] In the double polariza- 
tion case, either one initial and one final polarization are 
measured (such as e~ and p) or the polarizations of both 
final particles (p and p) are measured. In the first case 
the relevant observables (such as Q z q x o) require the ini- 
tial lepton to have longitudinal (±z) polarization, which 
is difficult to achieve experimentally. In the second case 
the initial e + e~ beams are unpolarized, and the longitu- 
dinal polarization of one final particle and the transverse 
polarization of the other must be measured. Determin- 
ing the p polarization may prove to be an experimental 
challenge. 

Of these two general possibilities, the most attractive 
"next experiment" beyond unpolarized e + e~ — > J/ip — > 
pp scattering may be a measurement of the differential 
cross section with unpolarized leptons and only the final 
p polarization detected. This will determine sinx, which 
specifies \ up to the usual trigonometric ambiguities. 

Another interesting experimental possibility is to re- 
solve the phase ambiguity in unpolarized e + e~ — > J/ip — > 
pp scattering through a study of the closely related re- 
action e + e~ — > J/ip — > AA, which has recently been 
observed by BABAR [l8| using the ISR technique. Since 
the ippp and V'AA couplings are identical in the SU(3) fla- 
vor symmetry limit, a determination of J /ip AA couplings 



The unpolarized angular distribution for the process 
e + e~ — > J/ip — > pp, measured recently by the BES 
Collaboration, is inconsistent with theoretical expecta- 
tions for a pure Dirac J /ippp coupling. In this paper 
we have derived the effect of an additional Pauli-type 
J /ippp coupling, and find that this can accommodate the 
observed angular distribution. The J /ippp Pauli coupling 
may significantly affect the cross section for the charmo- 
nium production reaction pp — ► ir°J/ip 7 which will be 
studied at PANDA. There is an ambiguity in determin- 
ing the relative Dirac and Pauli J/ijjpp couplings from 
the unpolarized e + e~ — > J/ip — > pp data; we noted that 
this ambiguity can be fully resolved through measure- 
ments of the polarized reaction. The most attractive po- 
larized process to study initially appears to be the case 
of unpolarized initial e + e~ beams, with only the final p 
(transversely) polarized. Alternatively, measurement of 
the required polarization observables may also be possi- 
ble using the time-reversed reaction pp — ► J/ip — > e + e~ . 
It may also be possible to use self-analyzing processes 
such as e + e~ — > J/ip — > AA to estimate the Dirac and 
Pauli couplings in the closely related J/ipAA vertex. 
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